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The Jordan frame action for general disformal theories is presented and studied for the first time, 
motivated by several unresolved mysteries that arise when working in the Einstein frame. We present 
the Friedmann equations and, specialising to exponential functions, study the late-time cosmology 
using dynamical systems methods and by finding approximate solutions. Our analysis reveals that 
either the disformal effects are irrelevant or the nniverse evolves towards a phantom phase where 
the equation of state of dark energy is —3. There is a marginal case where the asymptotic state 
of the universe depends on the model parameters and de Sitter solutions can be obtained. Our 
Hndings indicate that the metric singularity found using the Einstein frame construction corresponds 
phantom behaviour in the Jordan frame and we argue that this is the case for general disformal 
theories. 


I. INTRODUCTION 

The elusive nature of dark energy has prompted a the¬ 
oretical interest in the cosmological dynamics of scalar 
fields (see [T] for a review) as a mechanism for driving the 
acceleration of the cosmic expansion. With the exception 
of the simplest models such as quintessence [ai and 
k-essence [4], theories that include an additional scalar 
are alternative theories of gravity [5] (see [5] for a re¬ 
cent compendium of cosmologically relevant theories); 
they include additional degrees of freedom that couple 
to matter, resulting in additional gravitational strength 
(or larger) interactions. 

Many (but not all) modified gravity models can be 
written in the schematic form 

-S' = / d'^xv^ [i^g{9n>') + </>)], 

( 1 ) 

where the scalar field (j) is taken to be dimensionless. This 
action describes a theory of gravity in the so-called Ein¬ 
stein frame. Cg contains tensor self-interactions of 5^1/ 
through curvature tensors such as the Ricci scalar and 
L(f, contains scalar self-interactions. No direct couplings 
of the scalar to curvature tensors are present and the 
modifications of general relativity (GR) are encoded in 
the coupling to matter. In particular, test bodies do not 
move on geodesics of 5^^, the Einstein frame metric, but 
instead responds to the composite metric gfj.i^{gtj,iy,4 >)— 
the Jordan frame metric. It was shown by Bekenstein 
[ZlIH] that the most general theory of a scalar coupled to 
matter that preserves causality is 

gni^ = C{(j),X)gg_yPD{(t>,X)d^,(t>d,,(j)-, X = d^cfd ,,(}>. 

( 2 ) 

Indeed, it has been shown that theories where matter 
is coupled to metrics of this form are free of the Os- 
trogradski ghost instability [HUT!?]. C{(j),X) is known 
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as the conformal factor, and its consequences have been 
well-studied, at least when it depends on </> only. Con¬ 
sequently, D{(j>,X) has become known as the disformal 
factor and the term the disformal coupling to 

matter, or simply the disformal part of the metric. Any 
theory where D{(j),X) 7^ 0 falls into the class of disformal 
gravity theories. 

Disformal couplings are ubiquitous in fundamental 
physics. They arise in the low energy effective action 
of string theory [20] and are linked to galileons through 
probe branes moving in higher dimensional space-times 
[mill]. They also arise in the decoupling limit of massive 
gravity [23] . In the context of Horndeski theories mm, 
the most general scalar-tensor theories with second-order 
equations of motion, they are the most general transfor¬ 
mation that preserve the form of the scalar-tensor sector 
when C and D depend on (j) only m- These are the 
motivation behind a recent phenomenological study of 
disformal theories in several different contexts [23 HSl- 
Wf\ . with most attention focusing on the case where C 
and D depend on cf only. 

The parametrised post-Newtonian (PPN) parameters 
for this class of disformal theories were calculated in 
[551 155] . where it was shown that they are completely 
determined by the cosmological scalar cfo- For this rea¬ 
son, knowledge of the cosmology of these theories is vital 
for determining their viability. The first steps towards 
this were made by [54] , who used dynamical system tech¬ 
niques to classify the cosmological solution space in the 
Einstein frame with the goal of identifying models that 
passed solar system bounds. Several new fixed points 
were found but all had one problem in common: the Jor¬ 
dan frame metric becomes increasingly singular as the 
fixed point is approached, corresponding to the lapse in 
the Jordan frame approaching zero. This may represent 
several pathologies with the theory including the lack of 
a non-relativistic limit and a freezing of the proper time 
for observers, all of which we will discuss in more detail 
in the next section. 

Currently, it is unclear whether the singularity is debil¬ 
itating for the theory or if it is merely an artifact of work¬ 
ing in the Einstein frame. The two frames are equivalent 
for calculational purposes provided that the solutions are 
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interpreted appropriately when relating them to obser¬ 
vations. Furthermore, one cannot interpret observations 
until the proper time for an observer is aligned with the 
coordinate time and, since matter moves on geodesics of 
the Jordan frame metric, this frame is singled out for ob¬ 
servation^ [42]. Finally, Wetterich |43| has shown that 
frame transformations may introduce spurious solutions 
that solve the field equations in one frame but not the 
other. Given this, a study of the Jordan frame cosmol¬ 
ogy with a view to addressing these outstanding issues is 
certainly merited. 

This is the purpose of this paper. In the next sec¬ 
tion, we introduce disformal gravity theories in the Ein¬ 
stein frame and discuss the nature of the singularity, in¬ 
cluding the associated pathologies, in more detail. The 
applicability of the Einstein frame dynamical system to 
fundamental observers is also discussed. We next move 
on to study the Jordan frame cosmology. In section 
jlllj we present the Jordan frame Eriedmann and Klein- 
Gordon equations and use them to develop a dynamical 
systems approach to classifying the solutions for expo¬ 
nential scalar potentials and disformal factors. We find 
that theories where the disformal factor is small (in a 
manner to be made precise below) behave in a similar 
manner to quintessence but theories where the disfor¬ 
mal factor is large are not well described by a dynamical 
systems analysis in the sense that the fixed points reveal 
little about the late-time cosmology. Instead, we focus on 
finding exact solutions at late times in section [Vj Here, 
we show that the theory exhibits phantom behaviour at 
late times with an effective dark energy equation of state 
w = —3. 

Models that exhibit phantom behaviour are precisely 
those that suffer from singularities in the Einstein frame 
and thus we conclude that the singularity is indeed a 
physical pathology, the Jordan frame manifestation be¬ 
ing phantom behaviour. There is a marginal case that 
corresponds to a specific tuning in the parameter space of 
the theory. In this case, the asymptotic state of the uni¬ 
verse is a function of the model parameters and we show 
that it is possible to achieve asymptotically de Sitter solu¬ 
tions using a suitable tuning. We discuss our findings and 
conclude in section jvT] In particular, we argue that the 
qualitative features we observer here—quintessence fixed 
points and phantom behaviour—are features of general 
disformal models. 

Eor the reader interested purely in the cosmology of 
disformal models and not the singular nature of the dis¬ 
formal transformation, the Jordan frame cosmology is 
presented here for the first time and can be found in sec¬ 
tion |lll| onwards. The Eriedmann equations here can be 


^ Note that we do not claim that one frame is any more physical 
than the other, only that the Jordan frame is the frame where the 
coordinate time is aligned with the proper time for an observer, 
thus making any calculations simple to compare with other al¬ 
ternate gravity theories. 


used directly for computing quantities such as the lumi¬ 
nosity distance-redshift relation, which requires a trans¬ 
formation to coordinates appropriate for comoving ob¬ 
servers if one uses the Einstein frame. Furthermore, the 
non-phantom regions of the parameter parameter space 
of exponential models is presented in the conclusions (sec¬ 
tion VII where we also discuss the application to more 
general models. 


II. THE EINSTEIN FRAME 


In this section we present the Einstein frame action 
we will consider and use it to describe the singularity in 
the Jordan frame as well as the potential pathologies it 
presents. The action we will consider is 


S = j d^xv^Mpi^ 




The Jordan frame metric is 


BH4) 
Jv ■ 


9 ^. 1 / — 9 ^. 1 / T AO 


+ Sm[gtiiy] 

(3) 

(4) 


Specialising to the case of a flat Friedmann-Robertson- 
Walker (FRW) space-time: 

dsg = — dtg J- Ug da;g, (5) 


where we use subscript E’s and J’s to represent Ein¬ 
stein and Jordan frame quantities respectively, the Jor¬ 
dan frame singularity can be seen by computing the met¬ 
ric determinants [Uj: 


y=5 




_ / d</>y 

dtE j 


( 6 ) 


When E = 1 the Jordan frame metric is singular. Using 
Q, the Lapse in the Jordan frame is 

= I - S (7) 


so that 


dtj = N dtE- 


( 8 ) 


One can see that, cosmologically, the metric singularity 
corresponds to this becoming zero. 

There are several physical issues with the approach to 
the singularity. Eirst, the Jordan frame space-time is 

ds^ = —N'^ d^E + a(0^ (9) 


and so the proper time r for physical observers is |22| 


dr 

dtj 


= N, 


( 10 ) 


and so an observers proper time is frozen at the singu¬ 
larity. Furthermore, since = — I, the Lorentz 

factor is [42] 



( 11 ) 
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where v'^ = dxYdtj. Typically, one derives the New¬ 
tonian behaviour of the theory by looking at the limit 
n/c <C 1. In this case however, this is not sufficient. As 
demonstrated by [H] , one also requires S ^ 1 in order to 
have a sensible post-Newtonian expansion. Since the sin¬ 
gularity corresponds precisely to S —>■ 1, this behaviour 
is lost as the singularity is approached and there is no 
sensible Newtonian limit. The lack of such a limit was 
also noted by [33 |31j using Einstein frame coordinates. 
In this case, Newtonian quantities such as the total force 
diverge as E —>• 1. Since N is also the ratio of the speed 
of light to that of tensors, the Einstein frame interpreta¬ 
tion of this is that there are no particles that move with 
non-relativistic velocities in this limit. 

One obvious question is then: why not use FRW co¬ 
ordinates with unit lapse from the outset in the Jordan 
fame? In this case there is no apparent singularity at 
the level of the metric and any potential pathology must 
appear through the solution of the Friedmann equations. 
Indeed, our choice of coordinates such that is FRW 
is not a choice of space-time since no particles follow 
geodesics of Applying the change of time-coordinate 
([^ to Q one has 

ds^ = — dtj-I-a(tj)^ dx^, (12) 


where a(tj) = a(tj(tE))- This is an FRW space-time 
and so one can see that the singularity found taking 
to be FRW is simply a coordinate singularitj^ The 
one remaining question is that of a(tj). Currently, it is 
not known whether or not the transformation (|8| intro¬ 
duces any singularities into the spatial part of the metric. 
Said another way, is there some finite time tj such that 
a(tj) = 0? This is a difficult question if one begins in 
the Einstein frame. Equation Q is highly non-linear, 
and one requires an exact solution to provide an answer. 
Conversely, the Jordan frame is a perfect tool because 
one can classify the entire cosmological solution space 
using dynamical systems or other techniques. One then 
has the cosmological information that can be compared 
to data, as well as knowledge of any pathologies. One can 
identify the Jordan frame coordinates corresponding to 
the singularity found using the coordinates ([^ because 
applying the transformation ® one finds 


= 



1 - 

(N 

1+ 

vdoy _ 


-1 


(13) 


^ The reader should note that this is strictly true in the context 
of an isotropic and homogeneous cosmology. Whether or not 
there are other physical scenarios where a non-removable Jor¬ 
dan frame singularity is present is unknown, although we note 
that, to date, none have been observed. All of the pathologies 
that arise due to the singularity discussed in this section are the 
result of the cosmological singularity and the aim of this work 
is to understand the physical implications of this. For this rea¬ 
son, we focus entirely on the cosmological singularity and will 
not attempt to address the more general question of potential 
singularities elsewhere. 


and so in these coordinates the singularity corresponds to 
B d(/>/ dtj/A —)• oo. Another advantage of working in the 
Jordan frame exclusively is the following. Wetterich [35] 
has pointed out that spurious solutions can exist whereby 
a specific solution may be a solution of the Einstein frame 
equations of motion but not the Jordan frame equations. 
This potential problem is mitigated by working in the 
Jordan frame from the outset. 

All of the potential problems discussed above clearly 
motivate our study of the Jordan frame cosmology of 
disformal theories. Ultimately, we will see that when 
d(/)/ dtj —>■ oo i.e. in the limit where the singularity 
is present in Einstein frame time, the Universe under¬ 
goes strong phantom behaviour (w = — 3 for exponential 
models) and therefore the pathologies associated with the 
singularity are physical, the Jordan frame manifestation 
being precisely said phantom behaviour. In terms of the 
Einstein frame coordinates (i one can see that as the 
singularity is approached, the Jordan frame lapse tends 
increasingly towards zero and therefore, for comoving ob¬ 
servers, a large number of e-folds can pass in a small 
amount of proper time. When viewed in this manner, 
the phantom behaviour is hardly surprising. 

Before moving on to study the Jordan frame, we end 
this section by discussing the use of dynamical systems 
in both frames. The Einstein frame dynamical system 
for exponential models was studied by [HllMl- In order 
to achieve an autonomous system of equations, one uses 
the variable JVe = In ue as a time variable and chooses 
appropriate variables Xi that span the phase space of the 
system. Fixed points correspond to points in the phase 
space where 


dxj 

cINe 


= 0 


Vi. 


(14) 


It is assumed that Ie is a monotonic function of JVe and 
that Ne —>■ oo ^ Je —>■ cx), both of which are necessary 
for global attractors of the dynamical system to corre¬ 
spond to the asymptotic state of the system (see [33135] 
for the more technical aspects of dynamical systems the¬ 
ory). Note however that Je ^ oo does not necessarily 
imply that tj ^ oo. This depends on being able to inte¬ 
grate (|^ exactly and so fixed points in the Einstein frame 
do not necessarily correspond to the asymptotic future in 
the Jordan frame. One case where this can be achieved 
trivially is the case where V = 1 at the fixed point. In 
this case, disformal effects are absent and the theory be¬ 
haves in an identical manner to quintessence. Away from 
this limit, there are some important and physically rele¬ 
vant quantities that cannot be calculated using the Ein¬ 
stein frame. One pertinent example of this is Hubble 
constant 


H] = 


N ■ 


(15) 


In the Einstein frame. He ^ 0 at the fixed points and 
there are also fixed points where V —> 0 in the same limit. 
These are fixed points corresponding to the singularity. 
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This means that the behaviour oi Hj is undetermined. 
If the phase space were one-dimensional, one could sim¬ 
ply use L’hopital’s rule to find the asymptotic value but 
there is no higher-dimensional analogue of this theorem. 
For this reason, the asymptotic value depends on the 
phase space trajectory of the specific solution as it ap¬ 
proaches the fixed point. In this case, the Einstein frame 
dynamical system fails to achieve its goal of predicting 
the universal late-time behaviour since knowledge of the 
fixed points alone is not sufficient to know the asymptotic 
state of the Universe. This is an artifact of working in a 
coordinate system where the asymptotic state of dynam¬ 
ical system does not correspond to the limit of infinite 
proper time as seen by comoving observer^ Conversely, 
the Jordan frame dynamical system is perfectly able to 
predict all of the physically relevant quantities precisely 
because the coordinates are FRW with unit lapse from 
the outset. 

viewed as either a coordinate 


III. THE JORDAN FRAME 


From here on we work exclusively in the Jordan frame. 
For this reason, we will drop all unnecessary subscripts 
and tildes; it is to be understood that all quantities are 
Jordan frame quantities. The Jordan frame action is 
complicated compared with the simplicity of its Einstein 
frame counterpart, as are the derivation of the field equa¬ 
tions. For this reason, we give the calculation of the Jor¬ 
dan frame action and the field equations in appendix 
and present the final results here. 

We begin by defining the disformal coupling, 


d(j) 


(16) 


In general, /3 can be an arbitrary function of (f> but, in 
what follows, we set /? to be constant so that 


= 6^“^. (17) 

Second, the scalar potential is 

U(<^) = (18) 


where A is a constant and mg is an, a priori arbitrary, 
mass scale. These choices are made so that the equations 
exhibit a scaling symmetry that allows for the existence 
of scaling solutions and hence the dimension of the phase 
space is minimal [34]. With these newly defined con¬ 
stants, one can start building the disformal model. First, 


® One may wonder whether it is possible to choose Einstein frame 
coordinates to avoid this problem. Such a choice of coordinates 
would require working in a coordinate system where is part 
of the Einstein frame metric and would ultimately require one to 
mix frame variables in the equations of motion. 


we need the Friedmann equations: 
l2 

3H‘^ = Y + Vu + (19) 

=-y - 47rG'p„u^/^-b -b#) , (20) 

which are derived in Appendix |A 2[ The variable u is 
defined for notational convenience and is given by 


u = 1 -b 


5202 

A 2 ■ 


( 21 ) 


One can already see the advantage of working in the Jor¬ 
dan frame; the Friedmann constraint contains the disfor¬ 
mal scale A and so it is possible to compactify the phase 
space without using unphysical variables. This is in stark 
contrast to the Einstein frame, where the Friedmann con¬ 
straint is identical to that of the equivalent quintessence 
theory and it was necessary to use advanced techniques 
relating to fixed points at infinity to determine the late¬ 
time dynamics |34j . Next, we need the the scalar equa¬ 
tion of motion, which can be is expressed as: 




SnGpmB'^ 

A 2 


y3/2_|_yy2 _|_ 


/3S204 

A 2 


( 22 ) 

Because we are working in the Jordan frame, the scalar 
is minimally coupled to matter and one has the usual the 
continuity equation: 


aiJpm = 0. 


(23) 


Equations (19|-(23) contain all the necessary informa¬ 


tion about the dynamics of the system. The complexity 
of these equations make it impossible to find exact ana¬ 
lytic solutions and one method to analyse their late-time 
behaviour is to use a dynamical systems analysis. These 
methods are not new in cosmology. Indeed, they were 
applied previously to study quintessence m and disfor¬ 
mal theories in the Einstein frame [31]. Moreover, a dy¬ 
namical systems analysis is a powerful tool to calculate 
the late-time cosmological observables. The unfamiliar 
reader can find an introduction to dynamical systems and 
their use in cosmology in [I] [34| ]45l |46l l48H5f)] . 

Before proceeding to formulate the equations as a dy¬ 
namical system, we pause to discuss the physical observ¬ 
ables we wish to calculate using our subsequent analysis. 
Typically one is interested in the dark energy density 
parameter Ude and the equation of state w. These are 
difficult to define and several inequivalent effective vari¬ 
ables are often found in the literatur^B Instead of defin¬ 
ing effective quantities, we will look for quantities whose 


4 See |87l I514455| for discussions relating to effective quantities in 
scalar-tensor theories. 
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definition and interpretation are insensitive to the the¬ 
ory of gravitj0 One suitable quantity is the deceleration 
parameter 


which implies that 

;^=-(l + 9) (25) 

independent of the theory of gravity. In the case of 
icCDM, one has 

9=^(1 + 3weff), (26) 

where Wef^ depends on both Wm and wde- This moti¬ 
vates the definition 

= —1— (27) 

g < 0, or equivalently, w^b < 1/3 indicates that the 
cosmic expansion is accelerating and so we will use these 
classify the nature of the solutions. Formally, one may 
define 

IIde = 1 — flm, (28) 

and we will often refer to this quantity but the reader 
should be aware that this is not the same quantity that 
is inferred from cosmic microwave background or lumi¬ 
nosity distance measurement^ it is merely an indication 
of what is driving the evolution of the universe. 

IV. FORMULATION AS A DYNAMICAL 
SYSTEM 

In this section, we formulate the Friedmann-Klein- 
Gordon equations as a dynamical system and classify the 
fixed points. 

A. Construction of the Phase Space 

In order to make contact with the quintessence litera¬ 
ture, we begin by introducing the new variables 


^ By which we mean they describe properties of the FRW metric 
and are not found by comparing the Friedmann equations to 
those resulting from the Einstein-Hilbert action. 

^ We use the notation eff to denote a single composite quantity 
that describes the evolution of the universe and not an effective 
equation of state for dark energy. 

^ By this, we mean that the values of w and ^de ^re found by 
fitting the data to functional forms where w is constant, which 
is not necessarily the case for disformal models. 


Quintessence 



FIG. 1. The phase space of the system. The point (0,0,1) 
corresponds to matter dominated solutions and the edges of 
the cone correspond to dark energy dominated solutions. The 
phase space of quintessence coincides with the base of the cone 
located in the Z = 1 plane. 

where instead of differentiating with respect to coordi¬ 
nate time t, we differentiate with respect to V = In a(t). 
We denote derivatives with respect to N using a prime. 
This coordinate choice allows us to reduce the dimension 
of the phase space by one. Using these variables, we can 
rewrite the Friedmann constraint as 

I = + y'^u + ■ (30) 

As noted by [34], the disformal phase space is three- 
dimensional and so we require one more variable to close 
the system. [33] chose the variable z = BH/A but this 
results in a phase space that is non-compact. In particu¬ 
lar, note that u = l-|- 6 a;^z^ in this system so that neither 
the X- nor the z-directions are compact owing to the fact 
that Um can be arbitrarily small. Instead, one can work 
in a compact phase space by introducing the following 
variables: 

X = xu~i, Y = yu~^ Z = u~Y (31) 

When written in terms of these new variables, the Fried¬ 
mann constraint becomes: 

+ Y'^ + nm. (32) 

This implies that —1<A<1, 0<Y <1 and 0 < Z < 
1, therefore the phase space is compact. When written 
in terms of these variables, the phase space is a half-cone 
with the vertex located at (0,0, 0). This is shown in fig¬ 
ure Its base is the semicircle x"^ +y'^ + Um = 1 located 
in the plane Z = 1. This is precisely the phase space 
of quintessence and hence corresponds to the subset of 
the theory where disformal effects are absent. Any fixed 
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points that lie on the base of the cone therefore have 
late-time cosmologies that are identical to those found 
for pure quintessence theories with an exponential po¬ 
tential [IH]. Note, however, that their stability may be 
altered since the three-dimensional phase space implies 
the existence of a third eigenvalue, and that the other 
two eigenvalues may assume different values from those 
found in a purely two-dimensional phase space. Setting 
rim = 0 , one can see that the sides of the cone corre¬ 
spond to dark energy dominated solutions whereas set¬ 
ting X = Y = 0, one has flm = Z = 1, and so the point 
(0, 0,1) corresponds to a matter dominated solution. We 
therefore expect all physical trajectories to originate from 
its vicinity. 

Special attention must be paid to the tip of the cone 
X = V = Z = 0, which corresponds to what would 
be the metric singularity had we worked in the Einstein 
frame. This is a peculiar point because a fixed point here 
tells us absolutely nothing about the late-time cosmology. 
Typically, fixed points such as these indicate that the ef¬ 
fective dimension of the phase space is reduced and one 
typically requires centre manifold methods to find the re¬ 
duced phase space. Indeed, these methods are necessary 
for analysing the cosmology in the Einstein frame [34] . 
The reduced phase space is often unphysicaj^ and an 


alternate approach is to look for approximate late-time 
solutions given that one has some idea of which terms in 
the equations can be ignored at late-times. This is the 
approach that we will adopt in section [y] 

The case /3 = A/2 was identified bv |34j as a special 
parameter tuning in the Einstein frame where the dimen¬ 
sion of the phase space is reduced to two. This remains 
the case in the Jordan frame, where one has 

V{ct>)B\ct>)=ml (33) 


which implies a relation between X, Y, and Z\ 


y2 


A2 1 - Zt ■ 


(34) 


This relation is an additional constraint that must be 
satisfied and hence only two of the variables are inde¬ 
pendent. In terms of the three-dimensional phase space, 
the dynamics of the system are restricted to the two- 
dimensional surface where (341 is satisfied and hence the 
phase space is two-dimensional. Eor this reason the dy¬ 
namics of this case must be treated separately. 

Using equation (32) to eliminate U„i, equations (I9)~ 
( 22 ) can be expressed as a system of three autonomous 
first-order differential equations: 


dX 

dN 


dY 


dZ 

d/v 


y X'^(3- 9Zt^ -H 6X^Zi (y^ - 3Zt -k -k 3 (zt - 1^ (y2 


2Z2 (^X2 ^3^1 - -h (zt - l) (Z2 - y2)^ 

V6XZ -5zl^+ 2/3 (Zs - l) (U - Z){Y + Z)) - iV&PX^Z (zt - l) 


Y 


2Z2 (^X2 ^3^1 - l) -h (zt - 1^ (Z2 - y2)^ 

X^ (3 - 9Zt^ -h ex^zl (y^ - 3zi -h 2Z2) -k 3 (zt - l) z^f'^ 


2Z2 (^X2 (szi - 1 ) -k (zt - 1 ) (Z2 - y 2 )^ 

VQX^Z (a -f 6/3 (zt - 1 ) - 3AZt) - VgXX (zt - l) Z {2Y'^ + Z^) 


2Z2 (|x 2 ^' 3^1 - 1 ) -h (zt - l) (Z2 - y 2 )^ 
ax (zt - 1 ) (v^ ( 2 / 3 X 2 + Ay 2 ) _ 6 XZ) 

2Z4 (^x 2 (^1 - azi) -k (zt -(Z 2 - r 2 )^ 


Using equation (20), one finds 


(35) 


(36) 


(37) 


iU _ 

~H ~ 

X^ (dzi - 3 ) -f 4V6I3X^ (zt - l) Z - ex^zl (y^ - 2zi -k Z^) -h 2V6XXY'^ (zt - l) Z - 3 (zt - l) (r 2 _ z^)^ 

2Z2 (^x2 (1 - azt) -h (zs -1) (Z2 - r2)^ 

(38) 


In the sense that the variables are far removed from the underly¬ 


ing dynamical quantities such as H and <j). The phase space still 
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which can be used to calculate WeS and q. One also has 
^DE = ^ — in these variables. 


B. Fixed Points when fi ^ Xj2 


There are a total of five fixed points of equations 
(35)-(37) that we list in table [l] Table [IT] lists the 
interesting cosmological quantities at each point. The 
corresponding eigenvalues are listed below; only points 
(4) and (5) can be late-time attractors: 


(1) Cl — 62 — 63 — 0 

(2) 61 = 3, 62 = -2^6/3 - 6, 63 = 3 -b 

(3) 61 = 3, 62 = 2^6/? - 6, 63 = 3 - y^A 

( 4 ) 61 = A( 2/3 - A), 62 = i (A 2 - 6), 63 = A 2 - 3 

(5) e. = ¥-3, e3 = -j(l + ^a3V), 


e3 = -|(l- 


V24-7A2 

A 


Interestingly, not all solutions are fixed points. Point 
(1) is actually a fixed line, hence the zero eigenvalue. As 
discussed above, we will deal with this point using late¬ 
time solutions rather than dynamical systems. We note 
that X = Y = Z = 0isan independent fixed point not 
shown in the table. It corresponds to a matter dominated 
solutiorj^and so cosmologically viable trajectories should 
begin near this point. One can see that it is a saddle 
point and so trajectories will eventually leave its vicinity, 
signalling the onset of dark energy domination. 

Points (2) and (3) are unstable nodes or saddle points 
that correspond to non-accelerating phases and so we 
will pay no further attention to them. Points (4) and 
(5) are both located in the plane Z = 1, which, as dis¬ 
cussed above, corresponds to quintessence subset. These 
fixed points are hence identical to the points found if 
one considers quintessence with an exponential potential 
in GR. In particular, point (4) is the dark-energy dom¬ 
inated point that exists when A < ^/6. Point (5) exists 
when A > Vb and exhibits a matter-like behaviour with 
WeS = 0. Unlike the case of GR, these points are not 
always stable when they exist. Indeed, one can see that 
both are unstable when 2/3 > A. When this is the case, 
the only stable fixed point is at the tip of the cone and 
the dynamical systems analysis does not reveal anything 
interesting about the late-time dynamics. Examples of 
this are shown below in figures and If the theory 
was GR and quintessence, the models with A = I and 


contains all of the late-time trajectories. 

^ This is point (1) in |34| table I and the first bullet point in 
m, Appendix C when the conformal parameter a = 0 (in both 
cases). 


Name 

X 

Y 

z 

Existence 

(1) 

0 

0 

0 < Z < 1 

all 

(2) 

-1 

0 

1 

all 

(3) 

1 

0 

1 

all 

(4) 

A 

v/6 


1 

A < Ve 

(5) 

Vi 

A 

Vi 

A 

1 

any 


TABLE I. The fixed points of the system ( |35[ )-(|37[) when /3 yf 
A/2. 


Name 

H'/H 

q 

WeB 

^DE 

(1) 

3 

2 

1 

2 

0 

1 

(2) 

-3 

2 

1 

1 

(3) 

-3 

2 

1 

1 

(4) 

A2 

2 

1 - 2) 

1 (A^ - 3) 

1 

(5) 

3 

2 

1 

2 

0 

3/A^ 


TABLE 11. The cosm ological variables at the fixed points of 
the system (35h(37l when /3 yf A/2. 


A = 4 should approach fixed points (4) and (5) indepen¬ 
dently of the other parameters. These models are plotted 
in figure with toq = A = Hq and /3 < A /2 (the model 
parameters are indicated in the captions). One can see 
that these points are eventually reached after a brief ex¬ 
cursion into the domain Z < 1. In figure we plot the 
same models but instead choose P > A/2. Once can see 
that, in this case, both models now evolve towards the 
tip of the cone. 

One can then conclude that models with P < X/2 have 
late-time cosmologies that are identical to quintessence 
whereas those with /3 > A/2 exhibit drastically differ¬ 
ent behaviour. We will calculate this below in section 
[V|but we note here for completeness that points (2)-(5) 
are identical to those found by (34p^ The reason for 
this is that, as discussed in section]Tl[ these points all 
have Z = 1, which corresponds to Bp/A = 0 i.e. no dis- 
formal coupling. In this limit, the Einstein and Jordan 
frames are equivalent, and so are the coordinates used to 
parametrise the phase spaces. 


C. Fixed Points when /3 = A/2 


As remarked above, the phase space is two-dimensional 
when P = A/2. To see this, we can use equation (33) in 


These are points (2)—(5) in that reference, table I. 




































- 1.0 




FIG. 2. The phase space trajectories for models with A = 1 
(red) and A = 4 (blue). In each case {3 = 0.3, mo = A = Ho 
and the initial conditions are = 1, — 0. The 

initial value of Ni = Inoi and H{Ni) were chosen such that 
the universe begins in a matter dominated phase at redshift 
10 with Slni = 0.99999. 


the Friedmann constraint (191 to find 


1 — (Id" 2/r^)a;^ 3- y'^ 


1 + 2 ^ ^ 
yZ 


(39) 


with /i = mo/A. In this case, the phase space is elliptical. 
There are two routes by which one can proceed to analyse 
the fixed points of the system. The first is to reformulate 
the equations in terms of x and y and use the Friedmann 
equation to eliminate One can then find the fixed 
points of the two-dimensional system and proceed in the 
usual manner. The second is to continue to work in the 
three-dimensional framework and apply the constraint, 
which in our variables is (see equation (34|) 


= 


1 - Ft 


(40) 


Here, we adopt the second approach in order to make 
contact with the previous analysis. Substituting the con¬ 
straint (401 into equations (|^-(|^ in order to eliminate 


Y, one finds the equations in the reduced phase space, 
which are given in Appendix [B] due to their length. 

The resulting fixed points are given in table in with 
the corresponding cosmological parameters given in table 
JVj The eigenvalues are 

_ 3 

“ 2 ’ 


FIG. 3. The phase space trajectories for models with A = 1 
(red) and A = 4 (blue). In each case 13 = 3, mo = A = Ho 
and the initial conditions are those indicated in figure 


= 62 = 0, (42) 

= 62 = 0, (43) 

3 ^ v'(2 (A^ -b 18A2 - 72) y^ + 72- 21X^) 

4 4Av/(2(A2-6) h^ + 3) 

(44) 

= A2-3, e2 = -3+y. (45) 

One can see that the first three points are unstable and 
so we will ignore them from here on. The fourth point is 
a deformation of the stable spiral found when /3 < A/2 
and when A > -s/b- Its form is rather cumbersome but, 
by taking the limit A —>■ oo, one can see that the largest 
eigenvalue tends to zero from below and is therefore sta- 
blepl The fifth point is a deformation of the stable 
attractor found when /3 < A/2 and when A < -s/S. One 
can see that when y > \/2 and A > -s/S the only stable 
point is the tip of the cone. Just like the analysis of the 
case P ^ A/2, this implies that the trajectories approach 
a centre manifold at late times. Again, we will analyse 
this case by looking for an approximate late-time solu¬ 
tion. The three possible types of solution are shown in 
figure 1^ 


(2) 

6l 

(3) 

6l 

(4) 

e± 


(5) ei 


(1) 6i 


62 = 0 


(41) 


11 


Whether or not it is an attractor or a stable spiral depends on 
the values of /r and A. 
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Name 

X 

Y 

Z 

Existence 

(1) 

0 

0 

0 < Z < 1 

all 

(2) 

0 

0 

0 

all 

(3) 

0 

0 

1 

all 

(4) 



(l-2p3)3A 

p < ^, A > a/3 

(5) 

A a 3/4 

/^ A2 /2aA= 1 A"/" 

/2A=M= I A"/" 

\2 ^ 6 

V6 

V 3 6 ( A2-6 + 

( A2-6 + 3 j 

^ ^ 1+2^2 


TABLE III. The fixed points and lines when /3 = A/2. 


Name 

H'/H 

q 

Weft 

XIde 

(1) 

3 

2 

1 

2 

0 

1 - Z3 

(2) 

3 

2 

1 

2 

0 

1 

(3) 

3 

2 

1 

2 

0 

0 

(4) 

3 

2 

1 

2 

\2 

_1 

(5i-l)(l-2„»)-" + l 

(5) 

a2 

2 

\2 

2 ^ 

0 

1 


TABLE IV. The cosmological variables at the fixed points 
when P = A/2. 


V. LATE-TIME SOLUTIONS 

In this section we address models that evolve towards 
the tip of the cone by looking for approximate late-time 
solutions. These were the cases /3 > A/2 and j3 = A/2. 


A. Solution when /3 > A/2 


At late times, one expects that the field has rolled down 
the potential sufficiently such that cj) ^ 1 and flm 1. 
Writing the Friedmann constraint (19) as 


= 


1-h 


2J20-\)4> 


rriQe 


A2 


^ -I-rugC ^'^ + 8TrGpu^, 


(46) 

one can see that the final two terms are negligible com¬ 
pared with the term and so we have 




' A2 ■■ 


(47) 



FIG. 4. The phase space trajectories for possible solutions 
when P = A/2. The blue line tends to fixed point (4) and 
corresponds to a model with p, = 1, \ = 1. The red line tends 
to fixed point (5) and corresponds to a model with A = 1, 
p = 0.5. The black line corresponds to a model with p = 1, 
A = 10 and tends towards the tip of the cone. In each case 
A = Ho and mo was fixed using the value of p. The initial 
conditions are those indicated in figure 


This approximate solution is shown in figure and one 
can see that it matches very closely with the numerical 
solution. Next, we can make the same approximations to 
equation (20) to find 


Changing from coordinate time to A = In a we have 




A2 


= 3, 


(48) 


which is solved by 


= (49) 


2p-X 


2 mo 


H2 - 2 “A3® 




(50) 


which, when written using N as the time coordinate and 
applying the solution becomes. 


3^2 5/3-2 

2 ““^ A3® ~ 2N ■ 


(51) 
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Taking the logarithm of both sides, differentiating with 
respect to N and using the relation 


n ' 


-3-2- 


.U' 

ll 


(52) 


we find 


H' 2(/3 + A) 

H {2/3-X)N' 


(53) 


This is plotted in figure and one can again see that the 
approximation works very well at late times. 

We can see that when /3 > A/2 the universe will ulti¬ 
mately enter a phantom phase where H'/H tends to 3, 
although many e-folds must elapse before the asymptotic 
value is reached. That being said, it is not necessarily 
the case that a large change in N implies a large amount 
of coordinate time has elapsed. Indeed, recalling that 
H = dN/ dt, the coordinate time is 


t{N) = / 

JNi 


dN' 


(54) 


Since the lapse is unity, this is the proper time for co¬ 
moving observers. For non-phantom solutions such as 
the quintessence-like trajectories found in section |IVB[ 
H is a decreasing function of N and so t{N) is an ex¬ 
ponentially increasing function. The phantom solutions, 
on the other hand, have H{N) increasing exponentially 
and so t{N) is a slowly evolving function at large N. 
Physically, this means that one expects a large number 
of e-folds in a short amount of proper time, and so the 
asymptotic phantom state is reached very quickly. This 
behaviour is plotted in figure and can be understood by 
considering the Einstein frame. As the universe expands, 
the field begins to roll and disformal effects become in¬ 
creasingly important. If the field does not begin to slow, 
the Jordan frame lapse approaches zero and little coor¬ 
dinate time evolves, despite the fact that the scale factor 
and field are evolving rapidly. When viewed in this man¬ 
ner, phantom behaviour is a natural consequence of the 
disformal coupling. 


We end this section by noting that the solution (53) 
implies that q « —4, or, equivalently, WeS ~ —3. An 
equation of state this negative is in strong tension with 
observational data [551150] but it is not necessarily the 
case that this value is reached at the present time. In¬ 
deed, examination of figure reveals that, for the model 
studied there, the asymptotic value is not reached until 
far into the future. Whether a model predicts that the 
Universe is in the phantom phase at the present time, 
or that it will undergo one at some point in the future 
depends on the initial conditions and model parameters 
such as mg and A, which do not determine the asymp¬ 
totic state of the Universe but do control how quickly it 
is reached. For example, if one were to tune toq S> Hq 
the field will begin to roll early on in the Universe’s his¬ 
tory and one would expect phantom behaviour today. 
Conversely, tuning mg ^ Hq will result in the field being 


over-damped due to Hubble friction and the phantom be¬ 
haviour will only ensue far into the future. When fitting 
cosmological probes of the background expansion to data, 
Bayesian analysis will likely favour regions of parameter 
space where the phantom phase has not yet begun and 
so it is likely that cosmologically viable models can be 
found. Such an investigation would make an interesting 
topic for future work. 



FIG. 5. The evolution of (/)(A) found both numerically (red 
solid curve) and using the approximation ( |49| l (blue dashed 
curve). The parameters used were /3 = 1.4, A = 2 and mg = 
A = Ho- The initial conditions are those indicated in figure 

m 


H 



FIG. 6. The evolution of H'/H found both numerically (red 
solid curve) and using the approximation ( |49| l (blue dashed 
curve). The asymptotic value of 3 is shown using the black 
dotted line. The parameters used were j3 = 1.4, A = 2 and 
mo = A = Hq. The initial conditions are those indicated in 
figure [^ 
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/f„t(N) 



FIG. 7. The coordinate time as a function of N for 
both quintessence-like solutions (blue) and phantom solutions 
(red). The parameters used were /3 = 1.4, A = 2 (red) and 
/3 = 0.3, A = 2 (blue). In both cases mo = A = Hq. The 
initial conditions are those indicated in hgure[^ 


B. P = A/2 

When /3 = A/2 we can write the Friedmann equation 
as 

= {l + 2n‘^)^+ 8 ttGpu^. (55) 

Again, one expects that ^ 1 at late times but, un¬ 
like the previous, case, there are no factors of that 
become large in this limit. Instead, the second term is 
negligible and one has 


^2 = 


A (2/x2 + 1) (a + VA2 + 12^2 _ 6^ 


1 - 2 /i 2 


3(1 - 


(61) 


Note that one requires 2^^ > 1 in order for this type 
of solution to exist. If the converse is true the solution 


tends to fixed point (4) found in section IV C Given this 
constraint, one can see that solution (1) is incompatible 
with out assumption that (j) ^ 1 because (/>' = (5i < 0 and 
so only solution (2) is viable. A scaling relation such as 
this implies a definite prediction for the asymptotic state 
of the universe. Indeed, since is constant one has. 


using equation (52), 


H' 3A 

71 


(62) 


which implies 


9 = 


WeS = 


9A 


VA2 -b 12^2 _ 6 -p A 
6 A 

V'A2 -bl2Ai2 - 6 -b A 


— 4, and (63) 

- 3. (64) 


One can see that in this case the asymptotic state of 
the universe is a function of A and /i. Note that since 
(5i > 0, the universe cannot accelerate with H'/H > 3. A 
natural question is whether it is possible for the universe 
to achieve an asymptotic de Sitter state? Setting the left 
hand side of (62) equal to zero and using (61) one finds 


this is achieved when 


l«(l-|-2/i^)-and (56) 

6 

Ok+ (57) 


where the second equation comes from taking the limit 
1 in equation (20). Unlike the previous case, it is 
not possible to find an exact analytic solution but one can 
find late-time scaling solutions by looking for solutions of 
the form (j)' = 5i, = ^ 2 . Under these assumptions, 

one is led to two equations for 5i: 


3 — — (2/^^ 4“ l) ^1 4“ 3^2 

A (52 <^1 3(^2 

~ 2 2 2 ’ 


(58) 

(59) 


which have the solutions 


^ ^ \/A2 -I- 12^2 _ 5 _|_ y 

= - 

2 A (2Y -b 1) (a - a/A 2 -b 12^2 _ 0^ 

'^^“l-2/r2 3(1-2 /x2)2 

(60) 

( 2 ) =- 

a/A 2 -b 12^2 _ 6 -p A 


A = (65) 

As an example, we plot the evolution oi H'/H as a func¬ 
tion of N for the case /r = 1 (A = V2) in figure One 
can see that the universe does indeed tend to a de Sitter 
phase at late times. 



FIG. 8. H'/H as a function of N for a model with /r = 1, 
mo = Ho, A = y/2 and A = Ho- The initial conditions are 
those indicated in figure 
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VI. DISCUSSION AND CONCLUSIONS 

This paper has presented and studied the Jordan frame 
formulation of disformal gravity theories for the first 
time. The Einstein frame has been studied extensively 
and motivates this study for several reasons. First, there 
is an apparent metric singularity that previous studies 
have found, both numerically and analytically, to be ap¬ 
proached on cosmological scales when calculating using 
the Einstein frame formulation. This result has some 
pathological implications but, as discussed in section [T^ 
it is currently unknown whether or not is it a physical 
pathology or merely an artifact of working in the Ein¬ 
stein frame. This paper has taken the first steps towards 
answering this by studying the Jordan frame cosmology 
and looking for equivalent pathologies. Second, disformal 
transformations from the Einstein to Jordan frame do not 
preserve the lapse. This has the result that the proper 
time for observers in the Jordan frame is not aligned 
with the coordinate time, which makes the interpreta¬ 
tion of Einstein frame calculations difficult from a tech¬ 
nical point of view. The Jordan frame does not have this 
problem since the lapse is unity from the outset. 

The first part of the paper was dedicated to analysing 
the phase space of solutions using a dynamical systems 
analysis. We were successful in compactifying the three- 
dimensional phase space so that all solutions lie inside 
of the half-cone shown in figure [l] Interestingly, the 
phase space of the equivalent quintessence model (found 
by turning off the disformal couplings) coincides with the 
base of the cone, which allowed for transparent compar¬ 
isons with quintessence. In particular, any trajectory 
that terminates on the base of the cone has a late-time 
cosmology that is indistinguishable from quintessence, at 
least at the background level. The fixed points on the 
base of the cone correspond to those found in the Einstein 
frame by previous studies precisely because disformal ef¬ 
fects are absent and the time variables used to describe 
the dynamics in both frames are identical. Trajectories 
at the tip of the cone yield no information about the 
late-time cosmology and it was necessary to find approx¬ 
imate late-time solutions in order to discern the asymp¬ 
totic state of the universe. In this case, one can only 
relate the Einstein and Jordan frame time variables by 
integrating a non-linear relation, and it is here that the 
power of the Jordan frame formalism becomes apparent. 

The cosmological behaviour can be summarised con¬ 
cisely in the j3-\ plane shown in figure]^ When (3 < A/2, 
all of the fixed points lie in the quintessence plane and so 
the late-time fixed points are identical to those found by 
[3S] , although their stability is different due to the phase 
space being three- instead of two-dimensional. When 
/3 > A/2, the only stable fixed point lies at the tip of 
the cone and so it was necessary to look for approximate 
late-time solutions. These were found in section Ivl where 
we showed that the universe asymptotes to a phantom 
state where Weff = — 3 = 3) independent of the 

model parameters. One can see from the various figures 


that the pathological behaviour is typically reached in 
the future for universes that start from matter domina¬ 
tion and so it may be possible to reconcile the models 
with current observations. In particular, there are sev¬ 
eral model parameters, such as toq and A, that do not 
alter the position of the fixed points nor the stability. 
One would therefore expect a wide region in parameter 
space where the Universe is close to ACDM today but 
may undergo a phantom phase some time in the future. 
Such a model is not at odds with current observations. 
When fitting the model to cosmological probes of the 
background cosmology, it is likely that this region will 
be preferred by Bayesian fitting methods, although such 
analyses lie beyond the scope of this work. 

There is a marginal case given by /3 = A/2 where the 
phase space is reduced to two. In this case we found 
two late-time attracting fixed points that lie inside the 
cone and one that lies at the tip. By looking for late¬ 
time scaling solutions we derived the asymptotic value of 
H/H"^ for solutions that approach the tip and, in partic¬ 
ular, were able to show that by tuning the parameters, a 
late-time de Sitter phase can be reached. 

This solution deserves further comment in light of the 
cosmological constant problem. In order to achieve the 
solution required for an asymptotic de Sitter phase it 
was necessary to tune /3 = A/2 but this is not enough. 
One must further tune A and /r to values where the fixed 
point at the tip is the only stable one and the asymp¬ 
totic value of Weff is exactly —1. The theory does not 
contain any sort of protective symmetry and thus the 
tunings required are unlikely to be technically natural. 
Furthermore, the model has nothing to say about the 
old cosmological constant problem because we have set 
all contributions to the cosmological constant from both 
the scalar and matter sectors to zero from the outset. 
Given this, the asymptotically de Sitter cosmological so¬ 
lution found here has little to say about the cosmological 
constant problem, and the fine-tuned model is hardly a 
compelling alternative to ACDM. 

One of the goals of this paper is to discuss the metric 
singularity found by previous works using the Einstein 
frame formulation of the theory. The pathologies associ¬ 
ated with this singularity were discussed at length in sec- 
tionHU There, we noted that it is a coordinate singularity 
since one can find a gauge where the metric is perfectly 
regular and that it is apparently absent in the Jordan 
frame since one can work in this gauge from the outset. 
We showed that the singularity is located at the tip of 
the cone in this gauge and, furthermore, that trajecto¬ 
ries approaching the tip are those that exhibit late-time 
phantom behaviour. The physical manifestation of the 
singularity is then clear: the universe undergoes phan¬ 
tom behaviour a la m- Retrospectively, this is some¬ 
what to be expected from the Einstein frame behaviour: 
the approach to the singularity corresponds to the Jor¬ 
dan frame lapse approaching zero so that the clock for 
comoving observers slows down. A large number of e- 
folds can then pass in a short amount of time, which is 
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FIG. 9. The cosmological solutions found in this work. 


precisely the behaviour of a phantom Universe. 

We end by discussing the generality of our findings. In 
particular, the choice between a theory that is identical 
to quintessence, a phantom Universe or a finely-tuned de 
Sitter phase seems unappealing compared with simpler 
models. Here, we have only considered models where 
the scalar potential and disformal factor are exponential. 
This choice was made in order to yield the minimal phase 
space and preserve some of the scaling symmetry present 
in quintessence models. More general models will have a 
larger phase space that will require different variables to 
explore and one hence expects a new set of fixed points. 
Despite, this one would expect the qualitative features we 
have found here to apply. In particular, the fixed points 
were found to correspond to either phantom behaviour 
or the equivalent quintessence model except for a finely- 
tuned set of parameters. When written in terms of the 
cosmological variables and using N = In a as a proxy for 
time in the Einstein frame, the disformal coupling leaves 
the spatial component unchanged but the Jordan frame 
lapse is given by = 1 — B{(f>)‘^. All non¬ 
phantom Universes have —?► 0 at late times and so one 
expects a set of fixed points corresponding to the equiva¬ 
lent quintessence models precisely because the disformal 
coupling is set zero dynamically and the Jordan and Ein¬ 
stein frames are equivalent. One can then discern the 
requirements for the existence of new fixed points cor¬ 
responding to non-zero disformal couplings: either the 
Einstein frame Universe must be phantom so that H or 
(j)' increases without bound, or B{(j)) must be chosen such 
that it is a strongly increasing function of </>. This was the 
case with exponential models. Indeed, here we found that 
the disformal factor was only non zero for values of /3 that 
were large enough to compensate for the decreasing of H. 


One can then see that phantom behaviour is expected 
for any model where B{4') can increase rapidly enough 
at late times, what is not universal is the prediction that 
Weff = —3, which is likely to be a theory-dependent pre¬ 
diction. Said another way, one can design models that 
do not exhibit phantom behaviour by construction. One 
simple example of this is simply i? = I, which shows only 
quintessence fixed points. A more general example is the 
case of monomial potentials V{4)) ~ (^”, ~ 0 "* 

(with n and m positive even integers). In this case, one 
would expect (j) to roll to the minimum of the potential 
located at </> = 0 at late times so that B{(l)) tends to zero 
and the system behaves like quintessence. In light of the 
discussion above, we conclude that the general features 
found here—quintessence-like fixed points and phantom 
behaviour—are properties of more general disformal dark 
energy models. 

In order to find fixed points that were neither 
quintessence models nor phantom Universes it was nec¬ 
essary to fine-tune several model parameters to specific 
values. This corresponded to reducing the dimension of 
the phase space so that BH/A was fixed by the kinetic 
and potential energy. In this case it could neither grow 
without bound nor become zero. Such a fine-tuning is a 
very special property of the model considered here and it 
is unlikely to be a feature of more general models. More 
technically, the symmetries of the equations of motion 
~ and ~ ^{4’)) were crucial in 

allowing one to have the minimal possible phase space di¬ 
mension and to identify the requisite parameter tunings. 
Constructing other theories that exhibit these features 
would require looking at the symmetries present when 
the disformal coupling is absent and choosing the func¬ 
tional form of B{(j)) appropriately such that the dimen¬ 
sion of the phase space can be preserved with suitable 
parameter tunings. It is then clear that the novel fixed 
points found in the marginal case are not general and 
require finely-tuned simple models to exist. 

We have not included a conformal factor in our anal¬ 
ysis and it is unlikely that this will have any mitigating 
effects for the pathologies. Indeed, a conformal factor 
was included in the Einstein frame analysis of |34j with 
the only effect being to move the location of the fixed 
points. Such factors are strongly constrained by solar 
system tests and so the change is expected to be mini¬ 
mal. Indeed, if a conformal factor is present then 

the Cassini constraint on the PPN parameter 7 con¬ 
strains dlnA/d()) < 10“^ |l2l [63], which is the factor 
that appears in the equations governing the cosmological 
dynamics. 

Finally, one can relax the universal coupling and couple 
to dark matter only, at the cost of introducing violations 
of the equivalence principle, which are poorly constrained 
in the dark sector. In this case, the relevant fixed points 
for observers are those found by |34| , although the metric 
governing the motion of dark matter is that of a phantom 
universe and so one expects a drastic suppression of late¬ 
time structure compared with GR. 
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(Al) 


Bettoni and Liberati m have shown that the Horndeski 
action [24]— the most general scalar-tensor theory with 
manifestly second-order field equations—is invariant un¬ 
der disformal transformations and furthermore that the 
Einstein frame exists only when terms quintic in the 
scalar are absent. For. this reason, we expect that the 
Jordan frame action takes the form 


S = J d^xy^ X) + Gsi^, X)acl) + Gi{cj>, X)R{g) 

+G4,x [(n^)" - + S^[g^,], (A2) 

where Gi are arbitrary functions, X = —g^''d^(l}d,y(j)/2, 
and □ = g^'^ XpX y. Our strategy is then to transform 
each term in (Al) into the Jordan frame by inverting 


W a nd th en performing manipulations to get it into the 
fmm (A2). To accomplish this, we follow the methods of 
[mutr We begin by inverting (|^ to find 


Appendix A: Transformation to the Jordan Frame 

In this appendix we transform the Einstein frame ac¬ 
tion to the Jordan frame and derive the Friedmann and 
Klein-Gordon equations. 


1. The Jordan Frame Action 

Our starting point is the Einstein frame action (|^, 
which we write as 


5 /iid = gg.v ^2 and 

(A3) 

u. ~uu , B\<P) 

9 y A 2 -1 1 2 B 2 ( 0 )X’ 

i + A 2 

(A4) 

where all contractions are performed using g^^. 
introduce the tensor 

Next, we 

K^Ot _ -po; -pQ 

^ fllA ^ fll' 

(A5) 


A2(l+2iW) 

(A 6 ) 


using the identity 2 V[^Vi,]U^ = one finds [TS| 


S' = y d'^xMpi^ [Cg + Cff,] + Sm[ff] with 

£ and 

2 


= R%^.u + 2 V[^/c“ 


•AP ' 


O 


(A7) 


which, after making the appropriate contractions and us¬ 
ing equations ([^ and (A3), can be used to transform Cg: 






= \ n 


2B^{<i))X R 


B\ct>) 


2Ayi + ?W L 
BHcj,) 

A 2 ( 


(□</))" - 


A 2 2 
B{ct>)B^{cl)) 








+ 2 An(/> 


1 + 




A 2 






(A 8 ) 


This is not yet in Horndeski form; there are two quartic can remove this term and the quartic one by adding a 
and one cubic term that needs to be removed. Further¬ 
more, there is a term proportional to RgCX^ejXJ'^ (p. We 
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total derivative of the form with 

B\cj^) 


This simplifies the action to 




2 A Vl+^g^ 


^9 A , 2B^{4>)X R 

V + A 2 2 


. (A9) 


2A2,/i + ?^ L 


(□</>)" - V^V,</>V'^V"(^ 




A 2 ^/i^ 2 ^ L 


One can see that the quartic terms are in Horndeski derivative with 

form but there is still one cubic term that does not fit. 

This too can be removed by subtracting a second total 


C" = \ 1 


2B^(j))X B^{^) 
A2 B{(j)) 


2xa(j) + 

(AlO) 

(All) 


The action then becomes 


y=5 


= 1 1 + 


2B2(<^)X R 

A 2 ~2 


BH^) 


2A2 Jl + L 


(□ 0)2 - 


S^( 0 ) 


i?( 0 )A 2 / 


1 + 


2 B 2 ( 0 )x 

A 2 


1 + 


4S2(0)X 
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□0 


+ 2X 




4>4> 


BHcl>) 


2 B 2 ( 0 )a: 


1 + 


sl( 0 ) 


A2 


i?2(0)Jl+?iW 


(A12) 


which is now in the Horndeski form. 


Next, we need to transform C^j,. This is a lot simpler 
since one only needs to transform the metric determinant 
and the metric appearing in the kinetic term using (A4) 


to find 


7^9 


= t 1 + 


2 H 2 ( 0 )X 

A2 


X- 


2 h2 (0 )a:2 


A4 ( 


1 + 


2BR^)X\ 

A2 J 


The action is then in Horndeski form with 


B4.[^) 


G 3 ( 0 ,X) = -Mp 0 


B(«aAi + 'iSi^ 


1 + 


G2(0,X) =Mpi2Wl + 


2 H 2 ( 0 )X 

Is? 


2XB, 


4H2(0)X 
A 2 

‘^XBl{ct>) 
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2 H 2 ( 0 )X 2 


A4(i+2^) 


-y(0) 


-H(0) 

(A13) ■ 

(AM) 

(A15) 


(A16) 


Note that a similar action was obtained by [64]. 
2. The Field Equations 


reduce the action to minisuperspace using the coordi¬ 
nates 

ds^ =—N^{t) dt^ + a{t)^ dx^; 0 = 0(t). (A17) 

The Friedmann and Klein-Gordon equations can then be 
suiting field equations will be cumbersome and compli- found using the Euler-Lagrange equations for N(t\ a(t) 
cated. Since we are only interested in the homogeneous and 0(t) and setting N{t) = 1. Setting i3(0) = e^'^pjone 
and isotropic Friedmann equations it is simplest to first finds 


Given equations (A14|-(A16), it is clear that the re- 


S[N{t),a{t),m= [ dt— 

J /-i , 

V 4 + Af2A2 L 


o o , o p 


We do this for simplicity, it is not necessary to specialise at this 
stage but leaving the function general results in a far longer ex¬ 


pression. 
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-V{(t))N^ 



e2/30 N 

) 


2 6 ^^ ^ 

^ A2iV3 ^aN\ 7V2A2 ^ 


+ /37^+3 


A2oiV3 


(A18) 


where 


find equations (201 and (22). 


7 = 


1 + 


e‘^0'P 

2 _ 

iV 2 A 2 ' 


(A19) 


Appendix B: Dynamical System when /3 = A/2 


The Euler-Lagrange equation for N{t) yields the Fried¬ 
mann equation (after setting N = 1), which can be 
used in the Euler-Lagrange equations for a{t) and (j){t) to 


In this Appendix we present the dynamical system af¬ 
ter substituting the constraint (40) into (35)-(371 to elim¬ 
inate Y. They are: 


C7i = -A (4 + 3Z®/3 + 1 _ 4 /^ + 2V6XX^Z (2^? - 2 {2p? + l)zi + Z^/^ + l) 

-6X2^2 - 1 ) (^2^2 4^1 - 3 ) -h y/QXXZ^ (zt - - 3Z^ ~ 

t/ 2 -^ = AtA (3A^ (2fi^ + Zl -ij (2fi^ - 3Z5 + l) + 2y6A ( 2/^2 + 1 ) (zt - l) 

-fOA^Z^ (^Zt - 1 ) (2Ai^ -f 2Z3 - 3 ) - v^AAZ^ (^Zs - l) V SZ^ (zs “ > (B2) 

3A2 (v^AAZt - 6Z7/3 -t 6Z - v^A {2^^ + l) a) 

2 (a2 f-2/i2-A3Z3 - IJ - Zio/3-AZ2j 

where 

Ui = 2Z^ (zt - 1^ (x^ (2fi^ - 3Z3 YlJ +Z^ (zt - l)) (B4) 

U 2 = Z'^\j2-2Zt (zs - 1 ) (a^ (2y? - 3Z3 + + Z^ (z^ - l)) • (B5) 


Note that only two of these are independent since dif¬ 
ferentiating the constraint one has 


It is straightforward to verify that the dynamical system 
above indeed satisfies this relation. 


^X' _^Y' 4Z5 Z' 

A “ y + 3(1-Zt)^’ 


(B6) 
















